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FOREWORD 


The  work  reported  herein  conclude*  research  performed  at  the  Aeronutronic 
Division,  Ford  Aerospace  &  Communications  Corporation,  on  fluid-mechanic 
problems  involving  stability,  transition  and  turbulence  in  high-speed  shear 
layers  and  adjacent  streams  under  USAFOSR  Contract  F44620-75-C-0016.  During 
the  period  of  this  Contract  a  number  of  significant  experiments  were  done  on 
the  hydrodynamic  stability  of  hypersonic  boundary  layers,  resulting  in  a 
significant  fraction  of  the  existing  data  on  boundary- layer  stability  with 
and  without  wall  cooling.  Theoretical  and  experimental  work  was  also  done 
in  setting  up  rules  for  predicting  wake  behavior  in  gas-dynamic-laser  cavities. 
More  recently  the  effort  concentrated  in  formulating  a  general  transition 
theory  for  all  parallel  shear  flows.  The  latter  topic  is  the  subject  of  this 
report. 

The  work  reported  here  deals  with  the  formulation  of  the  theory  specifi¬ 
cally  for  free  shear  layers,  and  was  supported  by  the  Air  Force  weapons 
Laboratory,  Kirtland  AFB,  N.M.  (Dr.  P.  J.  Ortwerth) .  The  program  is  directed 
by  the  U.S.  Air  Force  Office  of  Scientific  Research  (Col.  Lowell  Ormand) , 
Bolling  AFB,  Washington,  D.C.  The  wind-tunnel  tests  were  performed  at  AEDC, 
Tallahoma,  Term.,  at  the  request  of  the  U.S.  Air  Force  Space  and  Missile 
Systems  Organization  (SAMSO/RSSE,  Maj.  M.  Sabin).  Mr.  J.  Donaldson  of 
ARO,  Inc.,  was  the  wind-tunnel  project  engineer.  Principal  Investigator  at 
Aeronutronic  was  Dr.  A.  Demetrlades  from  October  1978  until  May  1979. 

Dr.  A.  J.  Lade man  was  responsible  for  the  program  from  June  until  September 
1979,  during  which  time  he  contributed  material  for  inclusion  in  this  report. 

Part  I  of  this  report  presents  the  shear-layer  project,  while  Part  II 
presents  briefly  a  test  done  in  Aeronutronic' s  supersonic  wind-tunnel  on 
supersonic  turbulent  wakes.  The  list  of  symbols  and  Figure  captions,  as 
well  as  the  References  and  Figures  for  both  Parts  I  and  II  are  grouped 
together  in  the  beginning  and  the  end  of  this  report  respectively,  and  are 
numbered  sequentially  for  that  purpose. 
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ABSTRACT 


Part  1  of  this  report  presents  an  Investigation  of  the  transition  to 
turbulence  in  free  shear  layers.  A  theoretical  analysis  is  performed  to 
drive  a  formula  for  predicting  transition,  based  on  the  novel  approach  used  <t 

earlier  by  the  same  investigators  to  predict  transition  in  wakes  and 
boundary  layers.  The  formula  so  derived  shows  that  the  transition  Reynolds 
number  increases  as  the  non-dimensional  velocity  ratio  decreases,  and  also 
as  the  "fast  side"  Mach  number  and  "slow  side"  total  temperature  are  increased. 

An  experiment  to  test  the  theory  was  performed;  together  with  earlier 
transition  data,  the  new  data  provide  good  support  for  the  theory.  However, 
a  strong  dependence  of  the  transition  Reynolds  number  on  unit  Reynolds  num¬ 
ber  characterizes  the  available  data,  whereas  the  theory  does  not  predict 
such  dependence.  In  Part  II  of  this  report  a  brief  account  is  given  on  the 
measurement  of  temperature  fluctuations  in  a  cooled  compressible  turbulent 
wake.  The  results  endorse  the  novel  detection  technique  utilized,  and  support 
theoretical  predictions  made  earlier  under  this  contract. 
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1 . 1  INTRODUCTION' 


PART  I 

TRANSITION  IN  FREE  SHEAR  LAYERS 


The  author  of  this  report  has  devised  a  method  of  predicting  transition 
to  turbulence  in  shear  flows  such  as  wakes,  jet,  boundary  layers  and  free 
shear  layers.  The  method  is  based  on  simple  but  indisputable  physical 
principles,  and  thus  constitutes  a  necessary  condition  for  the  occurence 
of  transition.  Since  no  empiricism  is  involved,  there  are  no  "adjustable 
constants";  the  constants  appearing  are  well-defined  physical  properties, 
even  though  not  always  known  accurately  in  a  numerical  sense.  The  theory 
has  had  notable  success  in  wake  flows  (Reference  l)  as  well  as  in  boundary 
layers  (Reference  2) . 

Application  of  this  approach  to  free  shear  layers  (FSL)  is  needed, 
since  the  forecast  of  transition  in  an  FSL  is  the  first  step  in  mixing 
process  calculations.  The  practical  incentive  for  this  work  was  the 
necessity  to  predict  the  flow  in  chemical  laser  and  gas-dynamic  laser 
cavities,  where  such  mixing  is  generated.  There  is,  however,  a  more  impor¬ 
tant  objective:  a  general  method  of  transition  prediction  for  any  type  of 
shear  flow,  explicitly  including  effects  of  geometry,  compressibility,  etc. 
is  an  ideal  situation  because  it  can  be  applied  quickly  and  with  confidence. 
It  was  felt  that  the  past  success  of  this  technique  in  wakes  and  boundary 
layers  already  hinted  at  such  generality.  If  successful  in  shear  layers 
as  well,  the  method  could  indeed  represent  a  big  step  in  solving  the  long- 
vexing  problem  of  transition  to  turbulence. 

The  past  knowledge  on  FSL  transition  had  been  summarized  by  Birch  and 
Keyes  (Reference  3,  Figure  2).  Progress  seems  to  have  stopped  at  the 
definition  of  the  "transition  parameter"  Rex_  >.  and  the  empirical  finding 
that  this  parameter  increases  with  .  In  other  words,  all  we  knew  was 
that  the  transition  zone  (a)  scaled  as  Re  ,  (b)  moved  downstream  as  *'• 
decreased,  and  as  M,  increased.  Indicated  Improvements  were  the  formation 
of  a  rational  theory,  as  well  as  the  generation  of  more  experimental  data. 
Ihe  present  work  deals  with  these  two  issues. 

Beside  the  formation  of  the  theory,  the  present  work  supplies  new 
experimental  data  on  FSL  transition.  The  experiment  was  performed  in  the 
AEDC  Tunnel  B  at  hypersonic  speeds  although  the  levels  were  only  super¬ 
sonic.  The  so-called  velocity  ratio,  ,  was  verylsmall  (of  order  0.05) 
and  this,  together  with  the  Mach  numbers  achieved  (3  <  M,  <  4),  represented 
a  regime  previously  unexplored  in  FSL  transition  research.  The  data, 
joined  by  the  few  earlier  sets  of  data  available,  will  be  in  due  course 
compared  with  the  theory,  which  will  be  first  presented  below. 


1.2  TRANSITION  PREDICTIONS 


1 


1.2.1  METHOD 


The  theoretical  approach  to  predicting  transition  in  the  FSL  is  the 
same  as  utilized  earlier  by  Demetriades  (References  1  and  2)  for  deriving 
analogous  equations  for  the  wake  and  boundary  layer.  The  central  state¬ 
ment  of  this  method  is  that  the  turbulence  Reynolds  number 

i 

u  A 

ReA  ”  “vT  (D 

in  the  downstream  portion  of  the  turbulent  flow  following  the  transition 
tone,  has  to  have  a  minimum  value  Rey^  .  The  mechanics  of  prediction  then 
consist,  first,  of  assuming  that  the  particular  flow  under  scrutiny  is 
wholly  turbulent;  then,  Re^  is  computed  along  this  hypothetically  turbu¬ 
lent  flow.  The  point  along  the  flow  where  Re^  equals  Re^  is  the  transi¬ 
tion  point,  upstream  of  which  Re^  is  usually  smaller  than°Re^  .  Since 
no  turbulent  flow  can  exist  unless  Re^  >  Re^  ,  this  upstream  region  is 
laminar.  0 

Thus,  according  to  the  above,  the  present  task  consists  of  computing 
the  turbulence  Reynolds  number,  given  by  Equation  (1),  along  the  FSL.  This 
means  that  the  fluctuation  intensity  u',  scale  length  A  and  kinematic 
viscosity  u  should  be  computed  as  a  function  of  the  distance  x  from  the 
flow  origin,  going  downstream.  Before  this  computation  is  done  below, 
however,  it  is  necessary  to  clarify,  or  at  least  discuss,  the  following 
conceptual  difficulty. ,  The  quantity  Re^  varies  across  the  layer  as  well 
as  with  x,  and  since  u  is  zero  outside  the  FSL  and  has  a  maximum  in  it, 
then  Re^  will  also  have  a  maximum  in  it  according  to  its  definition 
(Equation  (1)).  The  question  now  is,  what  value  of  Re^  is  to  be  chosen 
at  each  x.  If  we  choose  the  maximum  Re^,  then  we  imply  that  the  flow  can 
locally  sustain  the  turbulence  if  Ke^  reaches  Re.  at  that  point  only. 

This  is  clearly  awkward:  assume,  for  example,  that  Re^  -  Re^  at  the 
center  of  the  flow  section;  turbulence  will  then  only  be  permitted  in 
the  exact  center  of  the  flow  and  not  anywhere  else.  Strictly  speaking 
this  is  an  inadmissible  claim. 

In  parallel  work  done  by  this  writer  to  investigate  transition  in 
a  boundary  layer,  the  same  question  could  not  be  circumvented,  and  the 
theory  had  to  consider  the  lateral  variation  of  Re^  at  each  x  station.  On 
the  other  hand,  in  applying  the  theory  to  wakes  (Reference  1) ,  successful 
predictions  were  made  using  the  Re^  magnitude  in  the  center  of  the  flow. 
Therefore,  despite  the  misgivings  of  the  preceding  paragraph,  we  will  here 
take  the  following  position:  if  a  turbulent  FSL  has  a  maximum  Re^  equal 
to  or  larger  than  the  threshold  Re^  ,  then  that  flow  is  permitted  to  be 
locally  turbulent.  If,  however,  conditions  are  such  that  Re^  at  its 
maximum  is  below  Re^  ,  turbulence  is  forbidden.  Transition  will  lie  at  the 
boundary  between  the  permitted  and  forbidden  portions  of  the  FSL. 

This  postulate  now  clears  the  way  for  computing  a  unique  Re^  (x)  for 
each  given  FSL.  We  know  that  u  will  have  a  maximum  on  the  dividing 
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streamline  (DSL);  thus,  all  quantities  In  Equation  (1)  must  be  computed 
there. 


1.2.2  CALCULATION  OF  THE  TURBULENCE  REYNOLDS  NUMBER 

The  objective  now  is  to  compute  the  variation  of  Re.  for  any  shear 
layer,  allowing  for  a  wide  variation  of  the  FSL  conditions  such  as  M. ,  l 
““d  ^02^01  ^n0  ®a8  differences  are  considered  here;  only  homo¬ 

geneous  flows  are  discussed).  In  line  with  previous  remarks,  only  the  Re. 
along  the  dividing  streamline  will  be  computed.  Thus,  for  each  FSL,  a 
unique  curve  Re^  (x)  will  be  obtained;  the  point  where  Re.  "Re.  will  be 
the  "transition  point".  Thus,  the  transition  distance  will  be  expressed 
in  terms  of  the  FSL  parameters: 


x_  -  x_  (M, ,  X ,  ) 

*  II  ’T01' 


(2) 


and  some  constants  which  will  be  discussed  in  due  course.  The  parameters 
forming  Re^  in  Equation  (l)  will  now  be  evaluated  one  by  one. 


1.2. 2.1  Kinematic  Viscosity 

The  subscript  convenction  will  be  to  use  "l"  for  the  faster  and  "2" 
for  the  slower  stream,  and  no  subscript  at  all  for  properties  on  the 
dividing  streamline  (or  "DSL").  Thus, 


u  -  (  —  )k  +  Kj 

v  T,  1 


where  k,  the  temperature-viscosity  exponent,  is  about  0.75  for  air.  The 
temperature  ratio  in  this  equation  can  then  be  found  in  the  Appendix: 

T  1  Y  -  1  2  ^02  ^  “  l)  /a\ 

t  "  2  < 1  +  1~rL  ><  i  +  > - r—  (4) 

1  1  T01  2(1  +  Y; 

1.2. 2. 2  Integral  Scale 

The  Integral  scale  A  is  perhaps  the  least  known  of  the  quantities 
needed;  in  contrast  with  the  Crocco  relation  giving  Equation  (4),  above, 
the  scale  comes  largely  from  making  an  educated  guess.  This  is  that  (a)  A 
is  constant  across  the  FSL  width  and  (b)  A  is  proportional  to  the  width: 


A  m  c  h 
1 
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The  width  h  is  known  from  the  work  of  Ortwerth  and  Shine  (Reference  4)  and 
others: 


h  -  g  (M  )  (6) 

o 

where  function  g(M.)  contains  the  "thickening"  effect  caused  by  the  Mach 
number,  and  0Q  is  the  incompressible  spreading  parameter  to  be  discussed 
farther  below. 

In  this  work  it  was  necessary  to  express  g(M_)  analytically  so  that 
transition  calculations  might  be  done.  Calculated  values  are  available 
from  Ortwerth  and  Shine  (Reference  4)  and  Oh  (Reference  5)  as  well  as 
experimental  data  from  a  number  of  sources  (see  Reference  5)  .  The 
differences  among  these  values  of  g(M,)  are  not  significant,  however,  and 
it  was  decided  to  compromise  by  curve-fitting  gO^)  into: 

g(ML)  -  0.3  +  0.7  exp  (-0.064*1^)  (7) 

1.2. 2. 3  Fluctuation  Intensity 

According  to  Ortwerth,  the  dividing-streamline  fluctuation  intensity 
is: 


u'  -  0.16  r^Xu,.  -  u2)  -  0.16  T(Ml)  ut 


(8) 


Like  the  function  g(M^) ,  ^(Mj)  was  for  convenience  approximated  by: 


f  (M^)  ■  exp  (-0.42Mj) 


(9) 


1.2. 2. 4  The  Turbulence  Reynolds  Number 

If  we  combine  the  above  equations  into  Equation  (l)  we  obtain: 
u,x 

)L  0.32  c,  : 

o 

The  critical  value  ReA  can  be  used  to  form  a  constant 


T,  k+l  x2 


ReA(x)  -  (  >£  0.32  ct  g-  ]  rOij)  gO^)  ( 


(10) 


Re* 


0.32  c.  =- 
*  o 


(ll) 
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Then  the  transition  Reynolds  number  is,  from  (10) : 


fcfl 


Re _ - 


XT  "  I'(M1)g(M1) 


(  “  ) 
Ll 


+  1 
7 r 


with 


Re 


XT 


XT  ■  transition  distance  from  flow  origins. 


(12) 


(13) 


This  is  the  desired  end  product  for  the  transition  distance.  It  depends  on 
X  not  only  by  the  factors  shown  on  Equation  (12),  but  also  implicitly 
through  the  ratio  T/T^  (see  Equation  (4) ) .  The  same  temperature  ratio 
also  includes  the  dependence  on  to2^Toi’  via  Equation  (4),  as  well  as  the 
implicit  dependence  on  M.  additional  to  the  factor  T(m  )g(M.)  in  Equation 
(12).  1  11 


The  relation  (12)  is  plotted  on  Figures  l  and  2.  In  these  graphs  we 
plot  Rex_/C  in  order  to  delay  discussion  of  the  constant  C  and  to  allow 
the  latter  to  adjust  itself  to  future  improvements  of  the  numerical  con¬ 
stants  composing  it.  The  plots  show  clearly  the  downstream  movement  of 
the  transition  point  (tone)  as  increases  and  as  X  decreases.  Qualita¬ 
tively  this  behavior  is  not  unexpected  and  agrees  with  earlier  and  current 
popular  notions  of  FSL  transition;  quantitatively,  however,  there  are 
substantial  differences,  as  will  be  discussed  below,  for  example  in  the 
dependence  of  transition  on  X . 


1.2.3  EVALUATION  OF  THE  NUMERICAL  CONSTANTS 

Accurate  knowledge  of  the  constituents  of  the  constant  C  of  Equation  (11) 
are  needed  for  numerical  application  of  the  transition  prediction  (Equation 
(12)).  In  previous  work  (see  Discussion  in  Reference  l)  the  threshold 
turbulence  Reynolds  number  Re^Q  had  been  found  to  be  approximately  15.  Use 
of  the  same  value  in  Equation  (11)  is  highly  desirable  because  it  will  test, 
in  the  long  run,  the  general  validity  of  the  present  approach  to  the  transi¬ 
tion  problem. 

No  information  seems  to  exist  for  shear  layers  on  the  constant  c. 
relating  the  integral  scale  and  the  layer  width  (in  Equation  (5)).  In  fact, 
what  is  needed  is  the  value  or  variation  of  c^  as  a  function  of  M. ,  X, 

T_»/Tqi»  etc.  including  its  change,  if  any,  across  the  layer.  This  serious 
snortcoming  will  be  met  here  by  assuming  c.  ■  0.2,  a  value  deriving  from 
wake  (Reference  6)  and  boundary  layer  studies  (Reference  l) . 

The  factor  0.32  in  Equation  (11)  derives  from  the  magnitude,  measured  by 
several  workers,  of  u'  at  its  maximum  point  in  the  shear  layer  (the  usual 
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finding  Is  0.16  (u,  • 
sequence  of  events;. 


u?) ,  and  here  a  factor  of  2  is  added  by  the  algebraic 
There  is  no  information  assuring  us  that  this  value  is 


unaffected  by  heat  transfer  (TQ2^'roi  ^  *  Actually,  experimental  data  on 

u  are  very  scarce  for  >  0,  as  well  as  for  sufficient  number  of  X  values; 
this  issue  is  certainly  far  from  settled. 


Using  the  most  commonly  known  values  of  c,,  Re^  ,  etc.  (the  spreading 
parameter  -  11.3)  we  can  then  compute;  ° 


C  - 


Re.  o 
Aq  o 


15  x  11.3 


0.32  n  c  0.32  x  3.14  x  0.2  - 


~  843 


(14) 


1.24  DISCUSSION  OF  THE  TRANSITION  EQUATION 

Equation  (12)  as  demonstrated  in  Figures  1  and  2,  and  augmented  by  the 
value  of  C  given  above,  presents  as  complete  a  set  of  predictions  as  ever 
attempted  up  to  the  present.  The  increase  of  ReXT  with  decreasing  X  is 
precipitous;  the  physics  of  the  problem,  of  course,  support  this  finding: 
at  X  -  0,  the  Re^  should  be  infinite  on  physical  arguments  alone.  When 
X  ■  1  (i.e.  u-  -0)  on  the  other  hand,  values  of  Rex  as  low  as  2C  (~  1700) 
are  obtained  In  the  common  case  of  incompressible,  adiabatic  flow  of  air 
(Figure  1) .  This  value  has  already  Increased  to  100C  by  the  time  X  has 
decreased  to  0.1  (Uj  *  0.82Uj). 

The  increase  in  Re^  is  also  very  fast  with  M  .  If  we  think  of  the 
case  u,  *  0  (X  ■  1)  Re^increases  by  a  factor  of  500  between  M.  ■  0  and 
Mj  -  5^1n  the  adiabatic  case.  It  should  be  stressed  that  this  is  the  first 
formal  enunciation  of  a  formula  on  FSL  transition  giving  quantitatively  the 
Mach  number  (compressibility)  dependence. 

This  is  also  the  first  occasion  in  which  the  transition  distance  is  shown 
to  depend  quantitatively  on  heat  transfer,  through  T  «/T  .  The  typical 

computation  of  Figure  2  is  done  for  three  values  of  TIT./T.:...  Shown  is 
the  result  that  as  the  stream  adjacent  to  the  "fast"  (subscript  "l") 
stream  is  heated,  the  FSL  transition  point  or  zone  moves  aft;  however, 
the  dependence  of  X?  on  T-./T..  is  not  apparently  as  strong  as  the 
dependence  on  M-  and  X .  A  noteworthy  feature  is  the  variation  of  Re—, 
with  X  when  T02'T01  “  1/3.  Here  8ee®s  that  for  certain  M^  values  the 
transition  sone  becomes  insensitive  to  X.  Computations  at  even  smaller 
values  of  Tq2^01  MOUl<*  worthwhile  in  this  context. 

There  is  very  little  previous  theoretical  work  for  comparison  with 
Equation  (12) .  Until  recently  it  was  thought  that  the  parameter  Re^X 
depended  only  on  M^,  i.e.  that 


Equation  (12)  shows  that  the  dependence  of  Re™  on  X  is  much  more  complex 
algebraically;  as  Figure  2  shows,  its  curves  nave  a  slope  varying  from  0 
to  1AZ  depending  on  and  X,  as  well  as  on  ^02^01 ' 


In  the  proposal  phase  of  this  work  (Reference  8) ,  this  author  had 
used  a  slightly  modified  version  of  the  present  approach  to  derive  a 
transition  formula  for  incompressible,  adiabatic  shear  layers: 


ReXT 


c 


X  +  1 


is  a  constant 


(15) 


Indeed,  for  TQ2/T ^  ■  1,  »  0,  Equation  (12)  reads 


ReXT  “ 


X  +  1 


(16) 


with  C  as  described  by  Equation  (11).  In  the  same  document,  the  author  had 
attempted  a  guess  at  the  extension  of  (15)  for  compressible  flows,  using 
the  then-available  test  data  of  Birch  and  Keyes  (Reference  3)  and  of 
Crawford  (Reference  9),.: 

exp  (1.52M.  +6.9) 

Rexr - jr1 -  <l7> 


This  expression  is,  thus,  essentially  the  locus  of  the  earlier  FSL  transi¬ 
tion  data  and  is  plotted  on  Figure  3  in  comparison  with  the  present  theory. 
The  agreement  is  surprisingly  good  and  acts  as  a  strong  endorsement  of 
Equation  (12),  as  well  as  it  assures  that  the  constant  C  and  its  ingredients 
(see  Section  1.2.3)  are  numerically  close  to  the  figures  quoted. 

Ortwerth  (Reference  4)  has  advanced  computations  of  Re™  which  follow 
the  spirit  of  Equation  (12) .  His  major  criterion  for  transition  is  closely 
related  to  the  critical  values  of  Re^  advanced  herein  (Reference  8)  and 
his  numerical  estimates  of  Re^  appea?  close  to  those  given  by  Equation  (12). 

A  comparison  of  Equation  (12)  with  the  data  obtained  in  the  present 
experiment  will  be  given  in  Section  1.5.7. 


1.2.5  QUALIFICATION  TO  THE  THEORY 

As  the  reader  is  well  aware,  the  theory  outlined  above  circumvents 
several  conceptual,  as  well  as  quantitative  issues.  A  partial  list  of 
these  conceptual  "short  cuts"  is  as  follows: 


(a)  Transition  is  supposed  to  occur  at  a  point;  actually,  it 
occurs  over  a  finite  region  which  may  be  quite  long. 


(b)  In  measuring  X^,,  the  virtual  origin  of  the  flow  Is  not  discussed. 
Generally  the  virtual  origin  is  not  known  "a  priori"  for  any 
given  flow. 

(c)  Transition  is  supposed  to  occur  when  the  Re;  on  the  dividing 
streamline  (or  the  maximum  Re^  at  any  flow  section)  reaches 
the  critical  value  Re*  .  This  is  equivalent  to  saying  that 
a  FSL  can  be  called  turbulent  even  if  turbulence  is  confined 
only  to  an  inf initensimally  thin  layer  along  the  dividing 
streamline,  a  clearly  unacceptable  concept.  The  application 
of  the  present  theory  to  the  boundary-layer  (Reference  2) 
discloses  that  the  consideration  of  the  off-axis  turbulence 
is  both  necessary  and  tedious. 

(d)  It  is  assumed  that  a  fully-developed  FSL  occurs  immediately 
after  transition.  Actually,  a  non-equilibrium  region  should 
exist  between  the  transition  tone  and  the  attainment  of  self- 
preservation. 

1 

(e)  The  theory  does  not  provide  for  a  dependence  of  7^  on  Re  , 
contrary  to  the  existing  observations. 

(f)  The  theory  disregards  initial  momentum  defects,  treating  the 
slipstream  trailing  a  shock  intersection  in  the  same  way  as 
two  boundary  layers  merging  downstream  of  a  partition. 

As  noticed,  the  theory  also  utilizes  information  drawn  from  the  turbulent 
state  of  the  FSL;  yet  this  information  is  woefully  inadequate  because  of 
the  lack  of  suitable  measurements.  Examples  are: 

(g)  The  theory  utilizes,  essentially,  Ortwerth's  theory  (Reference 
4)  on  the  variation  of  the  (maximum)  u  ,  which  is  in  turn 
verified  only  by  very  limited  data.  For  example,  u* (M. ) 

is  known  only  at  one  supersonic  Mach  number  (Reference1 10) . 

There  is  not  information  on  the  dependence  of  u  on  Tq2^toi 
or  on  \  at  supersonic  speeds. 

(h)  The  integral  scale  A  is  assumed  to  be  h/5  only  by  inference 
from  wake  and  boundary  layer  data.  There  are  no  appropriate 
measurements  of  A. 

(i)  The  kinematic  viscosity  is  seemingly  secured  by  the  use  of 
the  linear  Crocco  relation.  Yet,  both  in  boundary  layers 
(Reference  11)  and  wakes  (Reference  12)  the  linear  Crocco 
relation  is  now  known  to  represent  a  singularity  for  Prandtl 
number  equal  to  unity  only,  while  it  may  vary  considerably 

•  for,  say  TQ2/T01  4  l. 
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The  preceding  list  should  convince  the  reader  that  additional  work  is 
needed,  not  only  to  Improve  Equation  (12),  but  also  to  produce  the  experi¬ 
mental  data  on  which  that  formula  depends.  Recommendations  to  this  effect 
are  made  in  Section  1.7. 


1.3  EXPERIMENT  OBJECTIVE  AND  DESIGN 


The  objective  of  the  experimental  work  was  to  generate  a  free  shear 
layer  from  which  transition  data  could  be  obtained,  for  comparison  with 
the  theory  of  Section  1.2.  The  method  of  generating  the  FSL  became  an 
important  issue  at  the  outset  of  this  work.  In  practice,  most  FSL's  of 
Interest  are  those  between  two  co-flowing  streams,  e.g.  a  nozzle  dis¬ 
charging  a  stream  parallel  to,  and  into,  an  external  stream.  This  method 
was  shunned,  however,  because  of  the  interfering  effect  of  the  boundary- 
layer  growing  on  the  partition  separating  the  two  streams.  It  was  there¬ 
fore  decided  to  generate  a  "pure"  FSL  by  utilizing  the  slipstream 
emanating  from  the  intersection  of  two  oblique  shockwaves  in  supersonic 
flow,  a  method  previously  used  also  by  Birch  and  Keyes  (Reference  3) .  In 
contrast  to  the  experiment  done  by  the  latter,  it  was  decided  to  utilize 
two  planar  shocks  in  order  to  avoid  pressure  gradient  effects  generated 
when  one  or  both  of  the  shock  waves  is  curved. 


Considerable  work  went  into  the  design  of  the  method  for  intersecting 
two  shock  waves.  One  of  the  alternatives  studied,  for  example,  was  two 
opposing  wedges  resembling  a  supersonic  diffuser  inlet;  this  was  rejected 
in  the  end  because  the  "strength"  of  the  slipstream  produced  with  it 
(i.e.  the  difference  of  flow  speed  across  the  FSL)  could  not  be  raised 
without  choking  the  diffuser  inlet.  The  design  finally  chosen  consisted 
of  two  shock  waves  of  the  same  family  generated  by  a  double  wedge,  i.e.  a 
sharp-lipped  flat  plate  with  a  wedge  indentation  parallel  to  and  downstream 
from,  the  leading  edge.  The  flow  produced  by  this  arrangement  is  shown  on 
Figure  4.  Computations  were  made  to  see  what  strength  of  FSL  could  be 
produced  as  a  function  of  the  free  stream  Mach  number  and  the  angles 
and  o t  .  it  was  found  that  if  M,,,  <  6  or  so,  the  FSL  strength  was 
negligible  for  all  practical  angles  a.  and  ot  For  in  excess  of  8,  on 
the  other  hand,  the  FSL  strength  did  not  increase  much  beyond  what  could 
be  attained  at  ■  8.  Since  the  latter  was  also  the  design  Mach  num¬ 
ber  of  AEDC  Tunnel  B,  it  was  chosen  to  be  the  stream  Mach  number. 


Figure  5  shows  the  FSL  strength  attainable  at  Mach  8,  expressed  in 
terms  of  the  velocity  difference: 


as  a  function  of  and  a  Since  shock  separation  lies  beyond  the  bounds 
of  the  glot,  it  appears  that  the  desired  operation  would  occur  around,  say 
a.  ■  20  and  a  *  20  .  However,  wind-tunnel  choking  is  a  problem  when  the 
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frontal  area  of  the  model  placed  in  Tunnel  B  exceeds  about  200  in.  for 
the  available  indented-plate  model  (see  Section  1.4.2).  Furthermore, 
boundary- layer  separation  could  occur  if  «  was  large;  in  Figure  5  the 
chosen  and  should  stay  to  the  left  or  the  dotted  lines  which 
represent  estimates  of  maximum  allowable  a  before  separation  occurs. 

For  a  laminar  boundary  layer  on  the  base  plate  (Figure  4)  we  see  that 
Hankey's  criterion  (Reference  13)  precludes  all  but  the  smallest  ai  and 
angles.  For  turbulent  separation,  on  the  other  hand,  Reeve's  criterion 
(Reference  14)  implies  attached  flow  for  all  but  the  largest  or _  angles. 

This  criterion  and  the  "frontal  area"  curve  marked  "200  in.^"  enclose  the 
area  in  which  the  choices  of  and  or^  are  restricted. 

The  conditions  finally  chosen  are  shown  as  calculated  on  the  following 
Table  1.  Subscript  1,  as  per  the  usual  FSL  convention  of  marking  the  "fast" 
and  "slow"  side,  represents  the  fast  side  which,  according  to  Figure  4, 
lies  closest  to  the  plate.  According  to  this  Table,  if  the  indented  plate 
model  was  configured  according  to  the  last  two  lines  of  the  Table  at  M*  -  8, 
the  other  conditions  would  be  obtained,  regardless  of  the  PQ  and  TQ  of  the 
wind-tunnel.  However,  estimates  for  each  condition  shown  were  made  with 
Equation  (17) ,  to  ensure  that  the  available  PQ  range  would  cause  transition 
(for  each  condition)  to  appear  within  a  few  inches  from  the  origin  of  the 
FSL  (i.e.  the  oblique  shock  intersection) .  This  was  important  in  order  to 
utilize  the  main  diagnostic  tool,  which  was  the  tunnel  shadow  graph  system. 

The  experiment  design  also  sought  to  guarantee  the  attainment  of  the 
conditions  of  Table  I  by  insuring  two-dimensional  flow,  which  was  the 
basis  of  all  computations  made.  Specifically,  it  was  attempted  to  make 
the  model  as  wide  as  possible  in  order  to  avoid  edge  effects.  The  desired 
width  was  the  main  reason  for  making  the  model  large  and  introducing  the 
wind-tunnel  blockage  (choking)  considerations  referred  to  earlier.  A 
flat  plate  model  of  the  proper  width  was  actually  available  at  AEDC,  so 
that  model  fabrication  at  Aeronutronic  was  limited  to  modifications  needed 
to  add  a  ramp  (or  "wedge  plate")  to  this  model. 

The  experiment  design  was  implemented  at  AEDC,  which  provided  the 
wind-tunnel  and  test  support  services.  The  description  of  the  test  hard¬ 
ware  which  follows  is  drawn  primarily  from  AEDC  Report  TSR-79-V30 
(Reference  15) . 


1.4  EXPERIMENTAL  EQUIPMENT 
1.4.1  TEST  FACILITY 

VKF  Hind  Tunnel  B  is  a  closed-circuit  hypersonic  wind  tunnel  with 
a  50-in.  diam.  test  section.  Two  axisymmetric  contoured  nozzles  are 
available  to  provide  Mach  numbers  of  6  and  8  and  the  tunnel  may  be  operated 
continuously  over  a  range  of  pressure  levels  from  20  to  300  psia  at  Mach 
number  6,  and  50  to  900  psia  at  Mach  number  8,  with  air  supplied  by  the 
VKF  main  compressor  plant.  Stagnation  temperatures  sufficient  to  avoid 
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TABLE  I 


THREE  OPERATING  CONDITIONS  SELECTED  FOR  THE  TEST 


Condition 


air  liquefaction  in  the  test  section  (up  to  1350°R)  are  obtained  through 
the  use  of  a  natural  gas  fired  combustion  heater.  The  entire  tunnel 
(throat,  nozzle,  test  section,  and  diffuser)  is  cooled  by  integral, 
external  water  jackets.  The  tunnel  is  equipped  with  a  model  injection 
system,  which  allows  removal  of  the  model  from  the  test  section  while  the 
tunnel  remains  in  operation.  A  description  of  the  tunnel  may  be  found  in 
Reference  16.  In  this  experiment  the  tunnel  was  operated  at  Mach  8. 


1.4.2  MODEL 

This  model,  illustrated  in  Figure  6,  consisted  of  a  20- in. -wide  sharp 
flat  plate  with  a  ramp  attachment  positioned  8.5  in.  from  the  leading  edge. 
Brackets  located  between  the  ramp  attachment  and  the  flat  plate  were 
bolted  together  to  set  and  maintain  the  ramp  inclination  at  10- ,  15- ,  or 
20-deg  with  respect  to  the  flat  plate,  using  a  predetermined  bolt  hole 
pattern.  The  underside  of  the  ramp  leading  edge  was  bevelled  to  mate  to 
the  flat  plate  surface  when  the  inclination  angle  of  the  ramp  was  10  degrees. 
Spacers  were  provided  to  support  the  ramp  leading  edge  for  inclination 
angles  of  15-  and  20-deg,  No  seal  was  used  between  the  ramp  leading  edge 
and  the  flat  plate. 

To  minimize  flow  separation  in  the  compression  corner,  it  was  desirable 
to  maintain  a  turbulent  boundary  layer  ahead  of  the  flat  plate- ramp  junction. 
For  this  purpose,  distributed  roughness  elements  (No.  -40  grit)  was  used  as 
a  boundary- layer  tripped  mechanism.  Particles  were  applied  in  a  1.5  in. -wide 
band  parallel  to  and  0.5  in.  aft  of  the  leading  edge  of  the  flat  plate. 

Twenty-four  pressure  orifices  were  used  to  measure  model  surface 
pressures  and  monitor  any  lateral  pressure  gradients  existing  on  the  surface. 
The  locations  of  the  orifices  are  indicated  in  Figure  6. 

The  ramp  was  designed  and  fabricated  by  FACC  to  be  fitted  to  an 
existing  VKF  flat  plate  model. 


1.4.3  IN  STRUMENTATION 

Three  types  of  diagnostic  instrumentation  were  used: 

(a)  Pitot  tube,  for  flow  field  measurements  and  transition  detection. 

(b)  Shadowgraph  system  for  transition  detection. 

(c)  Model  surface  static  pressures  for  flow  field  measurements. 

Surveys  of  the  free  shear  layer  were  made  using  a  retractable  over¬ 
head  probe  drive  system  (X-Z  Survey  Mechanism)  designed  and  fabricated  by 
the  VKF.  The  mechanism  is  housed  in  an  air  lock  located  immediately  above 
a  port  in  the  top  of  the  Tunnel  B  test  section.  Access  to  the  test  section 
is  through  a  40-in. -long,  4-ln.-wide  opening  which  can  be  sealed  by  a 
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pneumatically-operated  door  when  the  mechanism  is  retracted.  Separate 
drive  motors  are  provided  to  (1)  insert  the  mechanism  into  the  test  section 
or  retract  it  into  the  housing,  (2)  position  the  mechanism  at  any  desired 
axial  station  over  a  range  of  35  in.  with  an  uncertainty  of  +0.01  in., 
and  (3)  survey  a  flow  field  of  approximately  10-in.  depth  with  an  uncer¬ 
tainty  of  +0.001  in. 

For  the  present  test  an  additional  drive  mechanism  was  attached  to  the 
foot  of  the  vertical-drive  (flow  field  survey)  strut  (see  Figure  7) .  This 
unit  was  a  rotary  actuator  with  the  probe  mounted  on  a  wheel  driven  by  a 
worm-gear  which  could  sweep  the  probe  tip  through  a  0.5-in. -length  arc  of 
5-in.  radius.  The  readout  potentiometer  of  the  rotary  actuator  was  cali¬ 
brated  for  the  full-scale  travel  (arc  length)  of  0.5-in.  and  was  equipped 
with  a  clutch  which  disengaged  the  poteniometer  drive  shaft  at  either  end 
of  the  travel.  The  clutch  permitted  the  probe  to  be  positioned  at  any 
desired  initial  location  without  affecting  the  resolution  of  the  readout. 
The  actuator  housing  afforded  storage  space  to  protect  probes  when  not  in 
use.  The  temperature  of  the  rotary  actuator  motor  was  monitored  through¬ 
out  the  test  to  insure  that  an  acceptable  operating  environment  was  main¬ 
tained. 

To  survey  the  relatively  thin  shear  layer  a  pitot  pressure  probe  of 
small  dimensions  was  required  to  avoid  size  effects  on  the  resolution  of 
the  profiles.  The  probe  was  constructed  using  32  mil  0D  and  22  mil  ID 
tubing  bent  to  facilitate  alignment  with  the  local  flow  direction  (see 
Figure  8) .  The  probe  was  flattened  at  the  tip  to  further  reduce  the 
probe  dimension  in  the  plane  of  survey.  The  resulting  lateral  dimension 
Increase  was  acceptable  since  the  flow  field  was  two-dimensional.  The 
probe  was  fabricated  by  the  VKF. 

Shadowgrams  of  the  flow  were  obtained  using  the  Tunnel  B  standard 
optical  system  operating  both  in  the  spark  (short  exposure)  and  continuous 
(time  exposure)  a lodes.  The  static  pressures  on  the  plate  surface  were 
measured  at  the  locations  Indicated  on  Figure  6. 


1.4.4  TL  MATRIX  AND  TEST  PROCEDURES 

The  text  matrix  shown  on  Table  II  indicates  that  the  experiment  design 
points  of  Table  I  were  fulfilled.  The  test  was  conducted  with  the  model 
placed  in  the  position  indicated  on  Figure  9.  Changes  in  the  angle  of 
attack  (l.e.  the  angle  <*, ,  of  Table  I)  we  accomplished  remotely  while  the 
tunnel  was  running;  to  change  (the  "ramp"  or  "wedge"  angle)  the  model 
was  withdrawn  from  the  flow,  to  enable  test  personnel  to  work  on  the  model. 

Test  data  in  the  form  of  shadowgraphs,  model  surface  pressure  measure¬ 
ments,  and  pitot  pressure  surveys  were  obtained  for  various  tunnel  con¬ 
ditions.  A  range  of  free-stream  unit  Reynolds  numbers  of  from  0.08-  to 
0.3*snlllion  per  in.  was  covered  by  varying  tunnel  pressure  P  between  200 
and  800  paia  while  maintaining  a  nominally  constant  tunnel  temperature  of 
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T  *  1320  R.  Increments  of  tunnel  pressure  of  50  psi  were  used  in  acquiring 
tRe  shadowgraphs  and  model  pressure  data  and  increments  of  100  psi  for  the 
survey  data. 

Shadowgraphs  of  the  shear  layer  were  obtained  for  each  test  condition. 
Standard  70-nsn  film  magazine  exposures  and  4-  by  5-in.  fast-developing 
exposures  were  made.  The  former  provided  a  file  of  the  shadowgraphs  and 
the  latter  facilitated  rapid  examination  of  the  optical  data  as  an  aid 
in  directing  the  course  of  the  investigation.  The  file  of  spark-shadow¬ 
graph  exposures  was  supplemented  for  certain  test  conditions  by  time 
exposures  made  using  the  continuous  light  source  of  the  optical  system. 

Model  surface  pressure  measurements  were  made  along  and  off  center- 
line  for  each  test  condition  to  assess  the  plate  finite-span  effects.  The 
ramp  surface  (centerline)  pressure  was  also  used  to  supplement  the  pitot 
pressure  measurements  obtained  from  the  surveys  of  the  shear  layer. 

Pitot  pressure  profiles  of  the  shear  layer  were  obtained  above  the 
ramp  near  the  trailing  edge  at  S  *  20  in. ,  approximately.  Each  survey  was 
begun  with  the  probe  immediately  below  the  shear  layer,  as  determined  by 
viewing  the  images  on  the  shadowgraph  screen  of  the  optical  system,  and 
the  initial  length  of  the  pitot  probe  (see  Figure  8)  was  aligned  parallel 
to  the  estimated  centerline  of  the  shear  layer.  Measurements  of  pitot 
pressure  were  made  at  25  to  35  positions  in  the  shear  layer  as  the  probe 
was  driven  in  discrete  steps  away  from  the  model.  The  locus  of  the  probe 
driven  by  the  rotary  actuator  was  an  arc  of  5-in.  radius  and  the  angle 
between  the  radius  and  the  centerline  of  the  shear  layer  was  typically 
40  deg.  The  readout  of  probe  position  was  limited  to  a  travel  of  0.5  in. 

An  arc  of  0.5-in.  length  with  a  radius  of  5.0  is  approximately  equal  to 
its  chord,  but  the  inclination  of  the  chord  with  respect  to  the  normal  to 
the  shear  layer  resulted  in  a  limitation  of  probe  travel,  with  readout, 
to  approximately  0.4  in.  in  the  normal  direction.  For  certain  test  con¬ 
ditions  the  "thickness"  of  the  shear  layer  at  the  survey  station  was  some¬ 
what  greater  than  0.4  in.;  however,  this  restriction  on  the  available  read¬ 
out  of  probe  position  was  not  considered  to  be  a  compromise  of  the  test 
objective. 

It  should  be  noted  that  the  roughness  elements  used  to  promote  boundary 
layer  transition  on  the  flat  plate  were  progressively  eroded  from  the  sur¬ 
face  during  the  testing;  however,  the  shadowgraph  pictures  of  the  flow 
field  do  not  indicate  any  resulting  separation  of  the  boundary  layer  at  the 
compression  comer. 


1.5  RESULTS 

1.5.1  FLOW  FIELD  MEASUREMENTS 

The  flow  field  achieved  for  each  configuration  was  important  to 
analyze  first,  because  it  shows  to  what  extent  the  design  flows  of  Table  I 
were  achieved. 
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TABLE  II.  TEST  MATRIX 


RAMP 

ALPHA 

PO 

TO 

RE (INF) 
-5 

GROUP 

ANGLE (DEG) 

(DEG) 

M(INF) 

(PSIA) 

(°R) 

x  10 

PER  IN. 

1 

10 

-15 

.0 

7. 

98 

397.6 

1333 

1.47 

2 

400.8 

1318 

1.51 

3 

399.4 

1319 

1.50 

4 

449.7 

1322 

1.69 

5 

500.4 

1320 

1.88 

6 

' 

f 

551.3 

1323 

2.06 

7 

7. 

99 

599.0 

1329 

2.22 

8 

1 

653.0 

1318 

2.45 

9 

7. 

98 

398.0 

1323 

1.49 

10 

4/ 

352.6 

1327 

1.31 

11 

7. 

96 

298.8 

1322 

1.13 

12 

7. 

95 

252.8 

1318 

0.962 

13 

7. 

94 

204.1 

1318 

0.779 

14 

8- 

00 

699.5 

1327 

2.59 

15 

752.3 

1313 

2.83 

16 

> 

t 

801.8 

1310 

3.03 

17 

-18 

.0 

800.1 

1322 

2.98 

18 

-20 

.0 

802.4 

1326 

2.97 

19 

747.5 

1325 

2.77 

20 

> 

r 

697.4 

1325 

2.59 

21 

7. 

99 

666.9 

1325 

2.48 

22 

653.5 

1322 

2.44 

23 

599.6 

1312 

2.25 

24 

\ 

1 

548.7 

1324 

2.05 

25 

7. 

97 

498.1 

1328 

1.86 

26 

446.0 

1322 

1.68 

27 

> 

f 

396.3 

1313 

1.51 

28 

-15 

.0 

399.2 

1326 

1.49 

29 

> 

r 

401.2 

1316 

1.52 

30 

7. 

94 

197.2 

1314 

0.756 

31 

7. 

96 

304.5 

1331 

1.14 

32 

7. 

97 

499.4 

1312 

1.90 

33 

7. 

99 

603.8 

1313 

2.28 

34* 

601.6 

1320 

2.25 

35 

> 

f 

' 

r 

600.5 

1312 

2.27 

36 

\ 

t 

8. 

00 

703.2 

1317 

2.63 

37 

15 

-15 

.1 

7. 

97 

401.1 

1321 

1.51 

38 

j 

L 

350.3 

1320 

1.32 

39 

-15 

\0 

J 

r 

303.2 

1321 

1.14 

40 

7. 

95 

247.3 

1325 

0.933 

41 

7. 

97 

449.9 

1329 

1.68 

42 

501.6 

1311 

1.91 

43 

7. 

98 

551.8 

1316 

2.08 

44 

7. 

99 

603.2 

1325 

2.25 

45 

I 

w 

649.3 

1320 

2.43 

46 

8. 

00 

701.4 

1324 

2.61 

47 

1 

752.5 

1321 

2.80 

48 

\ 

f 

\ 

1 

> 

V 

801.9 

1321 

2.99 

49 

10 

o 

CM 

1 

.0 

7. 

97 

401.1 

1323 

1.51 

50 

1 

501.2 

1323 

1.88 

51 

7. 

99 

601.2 

1317 

2.26 

52 

-19 

.9 

8. 

00 

703.1 

1330 

2.59 

53* 

804.2 

1323 

2.99 

54* 

801.3 

1322 

2.98 

55 

> 

1 

/ 

r 

801.1 

1323 

2.98 

*  Incomplete  Survey 


DATA 

TYPE 

Model 

Pressure 


i 

Model 

Pressure 

+ 

Flowfleld 

Survey 


y 

Model 
Pre  ssure 


y 

Model 

Pressure 

+ 

Flowfleld 

Survey 


Typical  plate  surface  pressure  measurements  are  shown  on  Table  III  and 
Figure  10.  Small  anomalies  in  the  pressure  distribution  are  insignificant 
and  thought  to  be  caused  by  operational  (e.g.  calibration)  errors  in  the 
data  acquisition  system.  It  is  seen  that  the  three  pressures  measured 
along  the  span  of  the  wedge  plate  are  very  nearly  equal,  indicating  little 
or  no  three-dimensionality  (edge  effects)  in  that  flow. 

The  measured  surface  pressures  are  for  the  three  test  configurations 
have  been  compared  with  theory.  Generally,  the  base  plate  pressures  were 
found  about  10%  higher  than  expected,  while  on  the  wedge  plate  they  lay 
below  the  expected  levels  by  an  equal  amount.  These  differences,  thought 
to  be  due  to  errors  in  the  angle  of  attack  setting,  are  at  first  glance 
appreciable;  but  when  the  angle-of-attack  errors  needed  to  produce  these 
differences  are  computed,  the  results  are  much  better.  It  is  estimated 
that  these  disparities  can  be  caused  by  misalignments  of  order  0.01  degree, 
i.e.  less  than  one  minute  of  arc,  which  is  certainly  within  the  error  band 
of  the  Tunnel  B  mechanism  at  AEDC.  It  is  concluded  that  the  observed  sur¬ 
face  pressures  meet  expectations  to  a  high  accuracy. 

The  pitot  pressure  surveys  across  the  FSL  combined  with  the  known 
tunnel  supply  pressure  and  the  static  pressure  (which  is,  for  every  data 
group,  identical  throughout  with  the  measured  wedge  or  ramp  surface  pres¬ 
sure)  makes  it  possible  to  compute  the  conditions  in  the  two  streams 
(regions  1  and  2  of  Figure  4)  bounding  the  FSL.  Such  conditions  include 
the  Mach  numbers  M^  and  M-,  velocity  ratio  u2^ui  etc-*  and  are  typically 
shown  in  Table  IV,  as  well  as  in  Figure  11.  ZIn  the  latter  it  is  seen  that 
the  design  Mach  numbers  M^  and  M^  are  best  attained  at  the  lower  tunnel 
pressure,  but  all  differences  are  anyway  so  small  that  one  can  safely 
claim  that  the  design  conditions  were  met.  Specifically,  Table  IV  states 
that  the  Mach  numbers  computed  by  the  theory  (in  the  design  stage)  were 
duplicated  within  l%-2%,  the  velocity  and  density  ratios  within  0.l%-0.3%, 
and  the  parameter  X  within  2.5%. 

It  will  be  noticed  that  Table  IV  does  not  include  the  case  ■  15°, 

■  15°  (Condition  1  of  Table  I).  The  flowfield  (and  the  transition) 
data  for  this  configuration,  represented  by  groups  37  through  48  in  the 
test  matrix  (Table  II)  were  anomalous  and  scattered.  This  problem  was 
explained  after  the  test  by  the  shadowgrams,  which  showed  that  the  FSL 
produced  by  this  pair  of  merged  and  therefore  interacted  with  the 

boundary  layer  on  the  wedge  plate.  This  condition  was  beyond  ghe  scope 
of  this  work,  and  therefore  no  further  analysis  of  the  *  15  ,  ■  15° 

case  was  made. 

The  discussion  has  so  far  dealt  with  the  experimental  verification  of 
the  edge  (boundary)  conditions  of  the  FSL;  of  course,  these  conditions  are 
very  important  in  the  FSL  transition  data,  to  be  presented  below.  The  flow 
inside  the  FSL  is  equally  important,  since  any  anomalies  in  it  would  cast 
suspicion  on  the  quality  of  the  flow  and  of  the  experiment  itself.  Good 
flow  quality  has  been  assured,  as  was  seen  by  inspection  of  the  documented 
data.  These  are  given  in  great  detail  in  the  Data  Package  appended  to 
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TABLE  IV. 


RESULTS  OF  FLOW  FIELD  SURVEYS 


Configuration  2 

(DEG)  15 

a 2  (DEG)  10 

(Average)  3.85 

(Theoretical  Expectation)  (3.92) 

M2  (Average)  2.95 

(Theoretical  Expectation)  (3«01) 

U2/U1  (Average)  0.922 

(Theoretical  Expectation)  (0.924) 

X.  (Average)  0.0404 

(Theoretical  Expectation)  (0.0395) 

P2 /p ^  (Average)  0.692 

(Theoretical  Expectation)  (0.691) 


3 


20 


10 


3.23 

(3.28) 


2.39 

(2.43) 


0.888 

(0.891) 


0.0591 

(0.0576) 

0.689 

(0.691) 
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Reference  15.  An  example  la  given  here  for  Group  35,  in  Table  III  and 
Figurea  12  through  15.  Note  that  beaide  the  Mach  number  diatribution  across 
the  FSL,  we  plot  the  ratio  M/M,,  u/u_,  T/T,,  etc. 


Figure  16  shows  an  Important  element  of  this  measurement,  which  is  the 
verification  of  the  "Crocco"  relation  across  the  FSL  (see  Appendix) .  The 
Crocco  relation,  given  in  Equation  (A. 1),  features  decisively  in  the 
transition  theory  of  Section  1.2;  if  it  were  Invalid,  the  results  of 
Section  1.2  would  be  Invalid  also.  Figure  16  shows  that  the  agreement 
between  the  data  (here  represented  at  random  by  Groups  35  and  51)  and 
Equation  (A.l)  is  indeed  very  good. 

In  summary,  the  flow  field  data  taken  with  the  pitot  probe  and  the 
model  surface  pressures  have  shown  that  the  desired  flow  was  achieved 
with  high  precision  for  the  two  configurations  2  and  3  of  Table  I .  This 
flow  can  next  be  used  to  interpret  the  transition  observations. 


1.5.2  TRANSITION  TO  TURBULENCE 


The  diagnostic  tool  for  turbulence  detection  was  the  shadowgraph 
system.  The  principle  of  detection  rests  on  the  fact  that  laminar  shear 
flows  usually  have  large  density  gradients  which  decrease,  due  to  enhanced 
spreading  and  diffusion,  when  the  flow  becomes  turbulent.  Since  optical 
systems  produce  a  response  proportional  to  derivatives  of  the  density 
they  can  be  used  for  transition  detection.  On  shadowgrams,  for  example, 
a  laminar  flow  would  appear  dark  (black) ,  and  transition  would  be  indicated 
in  the  region  where  the  flow  becomes  lighter.  Quite  frequently,  too,  the 
eddy  structure  of  the  turbulence  is  directly  visible  on  short-exposure 
photographs  taken  with  a  good  schlieren  or  shadowgraph  system.  A  trained 
observer  can  utilize  such  a  system  to  detect  transition  faster  and  more 
economically,  than  by  using  other  sensors. 


Figures  17  and  18  show  the  flow  field  as  observed  by  the  spark  and 
time-exposure  shadowgraph.  The  flow  field  follows  quite  faithfully  the 
theoretical  construction  of  Figure  4  and  the  zone  of  transition  in  the 
FSL  is  quite  clear.  The  time  exposure,  especially,  bears  out  fully  the 
hypothesis  that  the  FSL  "disappears"  after  transition.  The  transition 
zone  is  marked  on  each  picture  as  determined  by  visual  inspection;  note 
that  Figures  17  and  18  do  not  represent  the  same  conditions.  The  transition 
distance  X_  is  measured  from  the  FSL  origin  (the  crossing  of  the  two 
original  snocks)  and  the  transition  point  so  determined. 

AR0,  Inc.  developed  another  approach  in  processing  the  same  photo 
prints  from  the  shadowgraph  system,  which  consists  of  scanning  the  photo¬ 
graphs  such  as  shown  on  Figures  17  and  18,  with  an  optical  densitometer. 
Iso-density  contours  (here  "density"  refers  to  that  of  the  photo  print 
contrast  level)  can  be  plotted  directly  as  in  Figure  9,  page  22,  of 
Reference  15,  or  broken  up  into  profiles  of  optical  density  ("grey  level") 
such  as  shown  on  Figure  19.  Low  grey  levels  in  this  figure  Identify  regions 


seen  as  dark  or  black  by  the  viewer  of  Figures  such  as  17  or  18.  As  expected, 
the  dip,  which  in  Figure  19  represents  the  FSL,  gets  deeper  (darker)  as  one 
approaches  the  FSL  origin.  This  method  is  supported  by  automated  equipment 
at  ARO,  and  thus  recommends  Itself  for  this  transition  detection  problem. 

Only  a  limited  amount  of  work  was  done  with  this  technique  for  this 
test,  mainly  as  an  exploratory  evaluation  of  it  for  future  use.  Although 
based  on  the  same  principle  as  the  visual  inspection  of  photographs  (the 
recognition  of  varying  shades  of  grey),  the  densitometer  method  eliminates 
observer  bias  by  assigning  numerical  values  to  these  shades.  The  question 
therefore  is  whether  the  densitometer  method  agrees  with  the  visual  method 
employed  here  to  measure  X_  for  each  of  the  conditions  listed  on  Table  I. 

To  check  this,  the  densitometer  results  are  compared  with  the  visual  results 
on  Figure  20.  Transition  in  the  FSL  occurs  where  the  densitometer  trace 
departs  from  "black"  toward  "white".  It  is  seen  that,  for  this  group  of 
data,  this  occurs  exactly  where  the  visual  technique  also  put  transition. 

The  visual  method  is  thus  supported  by  the  densitometer  method. 


The  transition  distance  was  measured  visually  for  each  of  the  two 
spark  pictures  taken  at  each  of  the  conditions  tested  (i.e.  the  data  groups 
of  Table  II).  The  X^'s  listed  on  Table  V  and  plotted  on  Figure  21  show  an 
error  band  whose  two  extremities  are  fixed  by  the  two  photographs  at  each 
data  group;  that  is,  each  of  the  two  photos  gives,  generally,  a  slightly 
different  X_.  This  is  to  be  expected,  since  for  each  condition  transition 
usually  "jitters",  i.e.  oscillates  randomly  about  its  mean  position.  How¬ 
ever,  the  transition  Reynolds  numbers  on  Table  V  have  been  averaged  between 
the  two  data  points  for  each  data  Group. 


It  is  immediately  seen  on  Figure  21  that  X_  increases  (in  fact,' almost 
linearly)  with  P  .  This  is  surprising,  since  Section  1.2  indicates  a  con¬ 
stant  ReXT  regard less  of  P  .  Thus,  according  to  Equation  (12)  the  5L 
should  decrease  as  P  incriases,  an  expectation  common  to  all  shear  flows. 
Figure  21  says  that  Re  increases  rapidly  with  P  ,  both  because  P 
Increases  and  because  XI  increases.  According  to°Table  V  Re„_,  far  from 
being  constant,  varies  from  1  to  5  million  for  a  -  15°,  a  »I0o,  (\  -  0.04) 


as  P  varies  from  200Qto  800  psla;  Re„T  varies 
for  S  m  20  ,  ot  m  io°  over  the  same  P  range. 


ram  about  2  to  4  million 


The  unexpected  fact  that  Re^-  thus  varies  with  unit  Reynolds  number 
is  reminiscent  of  the  unit  Reynolds  number  effect  on  boundary  layer  transi¬ 
tion.  In  that  instance,  increasing  the  wind-tunnel  P  generates  increased 
acoustic  boundary- layer  radiation  from  the  tunnel  sidewall  and  accelerates 
transition,  spoiling  the  "constant  Re  ’’  concept.  Here,  however,  the  effect 
is  going  the  "wrong  way",  since  an  increase  in  P  increases  Re~.  The 
result  of  Figure  21  cannot  be  attributed,  theref8re,  to  the  customary  side- 
wall  radiation  effect.  The  series  of  shadowgrams  shown  on  Figure  22 
indicate  just  how  increases  as  P  increases. 
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TABLE  V.  MEASURED  TRANSITION 
DISTANCES 


X(INCH) 

FIRST  PHOTO  SECOND  PHOTO 


PO(PSIA) 


1.5.3  COMPARISON  OP  THEORY  WITH  EXPERIMENT 

The  first  concern  in  examining  the  present  data  centers  around  the 
unit  Reynolds  number  effect  just  discussed.  To  this  end,  we  present 
Figure  23  which  plots  the  present  Re  measurements  as  a  function  of  the 
unit  Reynolds  number  (equivalent  to  p  ) : 


As  already  discussed,  the  data  are  very  sensitive  to  Re  .  Plotted  for  com¬ 
parison  are  also  Crawford's  data  (Reference  9).  We  see  that  these  measure¬ 
ments,  too,  show  a  clear  dependence  on  Re',  and  that  they  furthermore 
increase  as  Re'  increases.  This  increase  is  not  as  pronounced  as  in  the 
present  test,  and  Crawford  actually  found  that  his  X~,  decreased  as  Re* 
increased.  Nevertheless  Crawford,  too,  has  apparently  found  an  Re1 
dependence  and  in  the  same  sense  as  the  present  measurements. 

The  theory  of  Section  1.2  does  not  reveal  a  dependence  of  Re  on  Re', 
and  hence  the  issue  remains  unresolved.  It  should  be  noted  that  any  theory 
in  which  X^  is  cast  in  terms  of  Re^  will  be  similarly  impotent  in  yielding 
a  unit  Reynolds  number  dependence.  xTo  avoid  speculation,  this  matter  will 
be  left  aside  and  the  present  data  will  have  to  be  discussed  further  as  an 
"error  bar"  in  the  Re^-,  M,  \  plane  rather  than  a  single  point.  Such  bars 
are  superimposed  on  tne  theoretical  predictions  of  Equation  (12)  in  Figure 
24.  The  experiment  is  lower  than  the  theory  by  a  factor  of  about  2;  how¬ 
ever,  the  absence  of  information  on  the  unit  Reynolds  number  effect  leaves 
this  comparison  unsettled.  For  instance,  if  P  in  the  present  test  could 
have  been  raised  by  a  factor  of  about  2,  complete  agreement  between  theory 
and  experiment  could  result. 

Figure  24  also  shows  the  data  of  Birch  and  Keyes  (Reference  3) ;  the 
agreement  with  the  theory  is  excellent.  Figure  25  compares  the  theory  (with 
1  ■  1  in  Equation  (12))  with  the  data  of  Crawford  (Reference  9)  .  The 
agreement  is  more  than  satisfactory  even  though  Crawford's  unit  Reynolds 
number  effect  casts  his  data,  too,  in  the  form  of  vertical  bars. 

It  is  obvious  from  the  above  that  the  theory  does  a  creditable  job  in 
predicting  the  vicinity  where  the  transition  zone  would  occur  in  a  shear 
layer.  Considering  the  range  of  data  (M^  from  2  to  6,  X.  from  0.04  to  1) 
the  comparisons  made  attest  to  a  large  improvement  which  the  present  theory 
causes  to  the  state  of  the  art.  We  should  stress,  of  course,  that  factors 
of  2  are  conmon  errors  in  transition  measurement  as  well  as  in  the  definition 
of  what  transition  is.  But  even  more  important,  the  qualifications  to  the 
theory  (Section  1.2.5)  should  be  kept  in  mind.  It  is  indeed  surprising 
that  the  theory  does  as  well  as  it  does  with  the  sparse  information  necessary 
for  its  formulation. 
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1.6  CONCUUSIONS 


1.  An  analytic  expresaion  has  been  derived  for  predicting  the  onset  o£ 
transition  in  homogeneous  free  shear  layers.  This  theory  is  based 
on  physical  principles,  utilizes  no  adjustable  constants,  follows 
the  author's  conceptual  approach  for  other  flows  (e.g.  wakes),  and 
divulges  the  transition  dependence  on  M.,  I  and  T.-/Tn  for  the 
first  time.  The  results  support  qualitatively  the  earlier  empirical 
ideas  of  the  dependence  on  M.  and  X ;  transition  is  found  to  move 
downstream  as  T02^T0l  lnc*eased' 

2.  An  experiment  was  done  to  enrich  the  data  base  with  which  the  above 
theory  can  be  compared.  The  test,  designed  to  produce  supersonic  shear 
layers  at  a  hitherto  unexplored  range  of  X  (X  <  0.1)  achieved  the 
designed  flow  conditions,  although  some  portions  of  the  matrix  could 
not  be  utilized. 

3.  The  agreement  between  the  new  test  data  and  the  theory  is  good;  the 
expected  large  downstream  movement  of  X_  with  greatly  decreased  X 
was  found,  although  there  is  still  a  numerical  disparity  by  a  factor 
of  about  2.  Earlier  data  from  two  other  sources  show  very  good 
agreement  with  the  theory. 

4.  A  strong  dependence  of  X^,  on  the  unit  Reynolds  number  has  been  found 
in  the  new  as  well  as  in  the  old  data,  which  cannot  be  rationalized, 
and  which  is  not  accounted  for  in  the  theory. 


1.7  RECOMMENDATIONS 

The  following  is  recommended,  based  on  the  discussion  of  the  theory  and 
on  the  conclusions  just  drawn: 

a.  Continue  the  development  of  the  transition  prediction  formula, 
Equation  (12),  with  the  purpose  of  explaining  the  observed 
dependence  on  unit  Reynolds  number. 

b.  Extend  the  theory  to  Inhomogeneous  mixing  layers  (different  gases). 

c.  Experimentally,  verify  the  observed  unit-Reynolds-number  dependence 
under  the  same  test  conditions,  but  utilizing  an  alternate  transi¬ 
tion  detector,  e.g.  the  hot-wire  anemometer. 

d.  Obtain  FSL  transition  data  for  regimes  not  probed  so  far  (e.g.  for 
0.1  <  X  <  0.2,  or  M1  >  5  or  both). 

e.  Measure  the  turbulence  intensity  u*  and  scale  length  A  for  any 
feasible  combination  of  X  and  M^  >  l. 


PART  II 

DENSITY  FLUCTUATIONS  IN  COOLED  SUPERSONIC  WAKES 

2.1  SUMMARY 

An  experiment  was  conducted  for  the  dual  purpose  of  measuring  the  in¬ 
fluence  of  nozzle-cusp  cooling  on  its  wake  in  a  simulated  gas-dynamic- laser 
cavity,  and  of  evaluating  a  new  technique  of  measuring  turbulent  fluctuations. 
The  turbulent  nozzle-cusp  wake  was  probed  with  a  dynamic  pitot  tube  which 
exhibited  the  wideband  response  expected  of  a  high-quality  turbulence  sensor. 
The  density  fluctuations  detected  as  a  function  of  cusp  temperature  agreed 
satisfactorily  with  theoretical  predictions. 

2.2  PURPOSE 

This  author  has  published  predictions  of  fluid  behavior  in  supersonic 
wakes  trailing  gas-dynamic-laser  nozzle  cusps  (References  17  and  18).  These 
predictions  address  the  laminar  wake  as  well  as  the  mean  and  fluctuating  flow 
in  cases  where  the  wake  is  turbulent.  Verification  of  this  theory  by  experi¬ 
ment  is  important  because  of  its  good  potential  as  a  design  tool  for  GDL 
cavities.  In  general,  the  agreement  between  this  theory  and  the  experimental 
data  on  the  mean  wake  flow  of  Peterson  and  this  author  -(References  19  and  2G) 
is  adequate  but  no  measurements  of  the  wake  turbulence  in  GDL  wakes  has  been 
reported  so  far.  Some  such  data  are  reported  here. 

2.3  EXPERIMENT  DESIGN  AND  TEST  SET-UP 

The  turbulent  wake  density  fluctuations  p'  were  chosen  for  measurement 
from  among  the  many  possibilities  partly  since  these  fluctuations  are  expli¬ 
citly  predicted  in  Reference  18  and  partly  because  of  their  practical  impor¬ 
tance  to  the  GDL  cavity  flow.  Specifically,  the  object  was  to  determine  the 
effect  of  cusp  cooling  on  p',  an  effect  which  the  theory  of  Reference  18 
treats  at  some  length.  The  test  set-up  was  identical  to  that  of  Reference  20, 
using  the  same  triple  nozzle  array  described  in  that  reference.  In  this  ar¬ 
rangement,  sketched  in  Figure  26,  three  continuous  streams  of  air  flow  each 
at  nominal  M  *  4,  generated  at  P  >*  755  mm  Hg  and  TQ  ■  606° R,  are  separated 
by  the  two-dimensional  wakes  she9  by  the  trailing  edges  of  the  two  nozzle 
cusps.  Detailed  quantitative  description  and  other  characteristics  of  this 
set-up  are  given  in  Reference  20  and  will  be  omitted  here  for  brevity. 

The  wake  examined  here  was  that  of  the  upper  cusp  (Figure  26)  which  can 
be  cooled  internally  by  circulating  L^.  Various  levels  of  the  cusp  wall 
temperature  T  can  be  achieved  by  throttling  the  coolant.  As  mentioned  in 
Reference  20,Wcooling  below  T  =  -50°F  is  possible  but  impractical  because 
tiny  ice  particles  then  form  on  the  throat  region  of  the  cusp  which  greatly 
distort  the  flow.  In  this  test,  measurements  of  p'  were  made  at  T  *  -40°, 

0°,  40°,  90°  and  126°F. 

■*  * 

Since  Reference  20  indicates  that  both  wakes  were  found  laminar,  or 
transitional  at  best,  steps  were  taken  to  ensure  turbulence  in  the  wake  of 
the  upper  cusp  by  tripping  its  boundary  layer.  This  was  done  by  completely 
covering  all  its  wetted  surfaces  with  No.  60  grit  deposited  on  a  thin  epoxy 
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layer.  Thla  scheme  was  successful  and  the  upper  wake  was  completely  turbulent, 
as  desired.  Thus,  the  effect  of  Tw  on  the  density  fluctuations  of  a  turbulent 
wake  could  be  studied. 


2.4  PRINCIPLE  OF  MEASUREMENT 

In  order  to  avoid  the  fragility  problems  and  protracted  data  reduction 
operations  with  the  hot  wire  anemometer,  this  measurement  of  p'  put  to  test  a 
novel  approach  to  turbulent  diagnostics,  which  Is  the  dynamic  pitot  probe. 

This  probe  has  been  used  elsewhere  in  the  past  only  in  a  qualitative  way  and 
without  the  benefit  of  a  rationale  for  its  use.  This  rationale,  given  by  this 
author  recently  (Reference  21),  is  as  follows.  The  instantaneous  fluctuation 
Ap  of  the  pitot  pressure  p  is  a  function  of  the  corresponding  fluctuations  Ap 
an!  Am  of  the  static  pressure  p  and  Mach  number  M,  respectively: 


API  “  AP  +  (7  _ 


7M2-1  M 


for  Y»  1.4;  here  Apr,  Ap  and  AM  stand  for  time  functions,  not  rms  or  ave¬ 
raged  quantities.  For  M  >2,  the  parenthesis  in  Equation  (18)  becomes  2, 
and  if  the  equation  of  state 

AP  Ap  *T 

—  mT+~  <19> 

is  also  used.  Equation  (18)  gives 

<20) 

»i  f>  " 

The  presence  of  the  velocity  f luctuations Au/u  in  (20)  can  now  be  elimi¬ 
nated  via  two  assumptions.  One  is  the  "standard"  assumption  of  hot  wire  ane- 
memetry: 

Ap  Ap  AT 

—  -O’— - x  (2D 

The  second  assumption  utilizes  the  author's  consistent  finding  (Refe¬ 
rences  22  and  23)  that  in  supersonic  wakes,  even  with  heat  transfer,  the  ve¬ 
locity  and  temperature  fluctuation  are  anti-correlated: 

AT 

—  -  -CY-l)  M2  ^  (22) 

^  u 

Inserting  (21)  and  (22)  into  Equation  (20) ,  we  find  the  gas  density  fluc¬ 
tuation  in  terms  of  the  pitot  pressure  fluctuations: 


I 


(23) 


where  primes,  indicating  rms  values,  have  now  replaced  the  instantaneous 
values  £p  ,  etc.  The  assumptions  underlying  Equation  (23),  mentioned  above, 
are  quite  safe  for  the  present  purposes  and  this  equation,  therefore,  formed 
the  basis  for  the  present  fluctuation  measurement. 


2.5  METHOD  OF  MEASUREMENT 


According  to  Equation  (23),  the  desired  goal  of  this  test,  i.e.,  the 
normalized  density  fluctuations  p'/p  ,  can  be  found  if  one  can  measure  the 
normalized  pitot  pressure  fluctuations  p'^/p^  and  the  local  Mach  number. 
Measurement  of  p'  ,  however  requires  a  sensor  of  proper  spatial  and  temporal 
resolution.  For  spatial  resolution,  a  sensor  size  on  the  order  of  57.  of  the 
wake  width  b  is  needed  or  about  0.1  cm  since  b  =  2  cm.  For  temporal  resolu¬ 
tion,  a  basic  rule  requires  a  frequency  response  extending  to  10u®  /b  or 
higher;  since  the  flow  velocity  is  on  the  order  of  7  X  10^  cm/sec  (Reference 
20)  the  response  should  extend  at  least  to  350  KHz. 

Until  recently,  pressure  sensors  meeting  these  requirements  were  diffi¬ 
cult  to  find.  A  relatively  modern  transducer  development,  the  Kulite  Model 
CQL-030-100,  was  chosen  for  the  present  work.  The  face  diameter  of  this 
sensor,  0.087  cm,  falls  within  the  size  requirement  quoted  above.  The  refe¬ 
rence  end  of  this  differential  transducer  was  hermetically  sealed  at  1  atmo¬ 
sphere,  and  the  sensor  was  repeatedly  calibrated  in  the  range  0  -  400  mm  Hg. 
abs.  This  is  a  small  fraction  of  the  transducer  100  psi  range  but  the  res¬ 
ponse  was,  nevertheless,  found  to  be  quite  linear.  Calibration  was  made 
easier  by  placing  the  transducer  permanently  with  its  tip  protruding  from  the 
front  of  a  cylindrical  holder  and  by  fitting  a  suction  cup  over  the  end  of 
the  holder  whenever  calibration  was  desired.  Since  frequent  calibrations  were 
needed,  this  scheme  made  them  possible  with  a  minimum  of  effort. 

Typical  calibrations  of  the  transducer  are  shown  in  Figure  27.  For  any 
single  calibration  the  output,  in  mv,  is  related  to  the  pressure  in  mn  Hg  by 


p  =»  A  +  Bp 
mv  mm 


(24) 


The  "zero  shift"  A  was  found  to  vary  from  one  calibration  to  the  next, 
making  this  transducer  unreliable  for  absolute  pressure  measurements.  The 
factor  B,  however,  remained  the  same  within  very  few  percent;  since  the  fluc¬ 
tuations  Ap  andAp  are  bv  the  differential  calculus  related  through 
mv  mm 
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it  follows  that  once  the  rms  voltage  output  fluctuations  p'  are  measured, 
the  fluctuations  in  pitot  pressure  can  be  deduced.  The  factor  B  was  found  to 
be  close  to  0.0067.  Together  with  Equation  (23),  Equation  (25)  completes  the 
formulas  needed  to  find  the  gas  density  fluctuations  from  the  measured  pitot 
tube  voltage  fluctuations. 

By  Equation  (23),  the  local  mean  Mach  number  is  also  needed,  which  was 
previously  found  by  measuring  separately  the  local  static  and  pitot  pressures 
(Reference  20).  For  expediency,  in  the  present  test  the  mean  pitot  pressure 
was  first  recorded  in  the  inviscid  flow  between  the  two  wakes  and  combined 
with  p  to  give  the  local  wake  edge  Mach  number,  and  thus  the  wake  edge 
static  pressure  p.  Using  the  finding  of  Reference  20  that  p  was  uniform 
across  the  wake,  it  was  next  assumed  that  this  p  was  the  same  on  the  wake 
centerplane.  The  measurement  of  p_  on  the  centerplane,  and  its  combination 
with  p,  thus  produced  M  *  M(0)  on  the  centerplane  for  insertion  in  Equation 
(23).  (Note  that  all  P'/P  data  taken  in  this  test  were  measured  on  the 
centerplane) .  In  the  beginning,  it  appeared  logical  to  use  the  same  Kulite 
transducer  to  take  these  two  values  of  p  ;  however,  this  was  impossible 
since,  as  mentioned  above,  the  erratic  values  found  for  A  in  Equation  (24) 
made  this  transducer  impractical  for  such  mean  (steady  state)  measurements. 
These  two  p  s  were,  therefore,  measured  with  an  ordinary  flattened  pitot 
tube.  1 


2.6  EXPERIMENT  GEOMETRY 

The  scope  of  this  test  was  limited  to  the  measurement  of  p'/p  on  the 
centerplane  of  the  turbulent  wake  for  the  few  temperature  ratios  Tw/T  men¬ 
tioned  in  Section  3.  The  dynamic  transducer  was  placed  at  x*5.8"  (15.7  cm) 

(Ax/hj-8)  downstream  of  the  trailing  edge  of  the  upper  nozzle  cusp  which, 
as  noted  earlier,  was  roughened  to  trip  its  boundary  layer  and  to  produce 
wake  turbulence.  The  flow  was  continuous  at  P0  =  756  mmHg  and  T0  *  606°R. 

The  probe  tip  lay  on  the  wake  centerplane  and  one  rms  voltage  output  for  each 
T  was  recorded.  After  this  was  done,  steady-state  (mean)  p-s  were  recorded 
sX.  the  same  conditions  by  an  ordinary  pitot  tube  which  was  also  used  to  mea¬ 
sure  wake-edge  pitot  pressures. 

2.7  RESULTS 

2.7.1  PERFORMANCE  OF  THE  DYNAMIC  PITOT  PROBE 

The  dynamic  pitot  tube  performed  very  well,  giving  a  large  a.c.  output 
in  the  turbulence.  The  key  criterion  of  its  operation  is  its  output  spectrum, 
which  is  shown  in  Figure  28.  As  noted  earlier,  signals  at  frequencies  of 
order  350  KHz  were  expected;  as  seen  from  the  spectra,  signals  at  frequencies 
as  high  as  640  KHz  are  detectable  with  little  trouble.  Furthermore,  pressure 
transducers  of  this  kind  have  a  "flat"  frequency  response  up  to  their  reso¬ 
nance  frequency;  one  advantage  of  the  CQL-030-100  is  that  its  resonance,  in 
consequence  :o  its  small  size,  lies  at  750  KHz.  Thus,  this  transducer  is  not 
only  amply  sensitive  to  the  p_  fluctuations,  but  has  a  flat  frequency  res¬ 
ponse  (an  advantage  over  the  not  wire)  over  the  entire  frequency  range  of 
interest  and  freedom  from  contamination  of  its  output  by  its  resonance  over 
that  range. 
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2.7.2  DENSITY  FLUCTUATIONS  IN  THE  TURBULENT  WAKE 

Table  VI  shows  the  results  of  the  quantitative  measurements  with  the 
Kulite  CQL-030-100.  For  nozzle  temperature  ratios  from  0.69  to  0.96,  the 
centerplane  Mach  number  changes  only  slightly;  its  effect  on  the  bracket  of 
Equation  (23),  however,  is  sufficient  to  convert  the  undecisive  variations  of 
Pp'  (from  8.5  to  9.3  nmHg  rms  in  Table  VI)  to  a  clear  variation  of  the  density 
fluctuations  p'/P- 

The  density  fluctuation  data  on  Table  IV  are  plotted  in  Figure  29  toge¬ 
ther  with  the  theoretical  prediction  taken  from  Reference  18.  According  to 
the  theory,  nozzle  cooling  suppresses  p'/p  by  about  277.  as  Tw  decreases  from 
Tq  to  0.6  T0,  and  by  about  19%  in  the  range  0.69  <  Tw/T  <  0.96.  In  the 
latter  range,  the  data  show  about  half  this  decrease  and  also  lie  about  307. 
below  the  theory. 

It  should  be  stressed,  of  course,  that  these  are  the  first  quantitative 
turbulence  data  reported  in  GDL  turbulent  wakes  and  the  first,  to  this  authorfe 
knowledge,  showing  the  cooling  effect  on  the  f luctuations.  At  this  stage, 
the  numerical  agreement  between  theory  and  experiment  should,  at  the  least, 
be  considered  satisfactory  and  the  evidence  of  turbulence  suppression  by 
nozzle  cusp  cooling  should  be  taken  as  conclusive. 

2.8  CONCLUSIONS 

1.  A  miniature  pressure  transducer,  employed  as  a  dynamic  pitot  tube,  has 
been  successfully  used  to  directly  measure  density  fluctuations  in  supersonic 
flow,  with  excellent  sensitivity  and  frequency  response.  This  performance, 
coupled  with  its  great  expediency  in  handling  and  signal  interpretation,  make 
this  technique  a  serious  rival  to  hot  wire  or  hot  film  anemometry. 

2.  It  has  been  verified  that  when  a  GDL  nozzle  cusp  is  cooled  below  the 
stagnation  temperature,  the  density  fluctuations  on  its  turbulent  wake  center- 
plane  decrease.  Qualitatively  and  quantitatively,  this  decrease  is  similar 

to  that  expected  by  the  theory  presented  by  this  author  in  an  earlier  publi¬ 
cation. 

2.9  RECOMMENDATIONS 

Recommendations  are  offered  here  both  as  regards  the  dynamic  pitot  tube 
technique,  and  also  as  regards  the  measurement  of  turbulence  in  cooled  wakes. 
The  quantitative  knowledge  of  turbulence  "suppression"  or  "amplification"  as 
the  heat  input  varies,  is  a  very  important  problem  at  this  juncture.  Since 
the  particular  wake  generators  used  could  not  be  cooled  below  a  certain 
limit,  it  is  recommended  that  the  test  set-up  is  changed  to  one  where  much 
lower  values  of  Tw/To  can  be  reached.  Specifically,  it  is  recommended  that 
the  test  is  repeated  in  the  same  wind-tunnel  at  Mach  3,  but  with  the  wake 
generator  consisting  of  a  thin  circular  tube.  Passing  any  amount  of  LN_ 
through  this  tube  should  be  quite  simple,  and  the  Tw/To  should  reach 
considerably  lower  values  than  in  the  test  discussed  here. 


Further  experimentation  with  the  dynamic  pitot  tube  is  strongly  encou¬ 
raged.  This  author  considers  this  so  Important  that  he  recommends  separate 
funding  of  a  program  to  measure  density  fluctuations  with  this  instrument. 
Specifically,  parallel  measurements  with  this  method  and  the  hot-wire  ane¬ 
mometer  should  be  done,  so  that  the  density  f luctuations  can  be  simultaneously 
measured  with  these  two  methods  for  purposes  of  comparison. 
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TABLE  VI 


FLUCTUATION  DATA 


Pq  ■  756  nmHg.  abs.. 

T  -  606° R,  x  - 
0 

5.8",  Ax/hj 

-  7.9,  y  - 

0 

CASE 

I 

2 

3 

4 

5 

T  (°F) 
w 

-40 

0 

40 

90 

122 

T  /T 
w  o 

0.693 

0.759 

0.825 

0.907 

0.96 

PT  (0)  (ran  Hg) 

58 

60.7 

58.9 

57.1 

51 

M  (0) 

2.185 

2.125 

2.195 

2.175 

2.10 

PT'  (0)  (ran  Hg) 

8.58 

9.36 

9.36 

9.36 

8.58 

P'  (0)/p  (0) 

0.0722 

0.0763 

0.078 

0.0797 

0.0789 
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FIGURE  4.  Schematic  of  the  model  and  flow  regime  nanenclature. 


FIGURE  7.  Schematic  o£  the  actuator  and  pitot  probe  placement. 


FIGURE  8.  Details  of  pitot  probe  tip. 


FIGURE  10.  Typical  pitot  transverse,  normalized  to  the  pitot  pressure 
in  region  2  (Group  35) . 


■■■■ ■■■■ ■■■■ ■■■■ ■■■■■■■■ ■■■■ 
■■■■■■■■■■■■■■■■■■■■■■■■■■a 
^^■^^hhmbubbbbbbbbbbbb?. 

liiiiiinilii 

■BBBBBBBBBBUflB 
■■■■■■■■■■■■■■ 

uiiiiiiiiiiiiniiiiSiiisiiii 

miiiiiiiiiiiiiiiiaiiiiiiii 

■■■■■■■■■■■■■■■■■■■■■■■■■■■a 

BBBBBBBflflflBBflflflBBBBBBBBBBBBB 

iiiignaniiiiiiiiiiiiJiiiii 

8S8S8S888S8SSSSS8SSSBS8888SS 


.it 


nasiiRBisai^iiiaiBRi 

■asEsassEginsssasHi 

!8SaE8E3EE8EE8S5BSaSSaSSa88! 


ess*  i 


h] 


FIGURE  14.  Variation  of  the  static  temperature  across  the  shear 
layer,  normalized  with  the  temperature  In  region  2 
(Group  35) . 


FIGURE  15.  Variation  of  the  density  across  the  shear  layer,  nor 
malized  with  the  density  in  region  2  (Group  35). 
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FIGURE  17.  Spark  shadowgram  of  the  flow. 
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FIGURE  20.  Comparison  of  transition  detection  between  the  denslto 
meter  and  the  visual  methods. 
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FIGURE  21.  Variation  of  the  transition  distance  with  tunnel  pressure. 
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FIGURE  22.  Shadowgrams  of  the  free  shear  layer  at  (from  top  to 
bottom)  200,  300,  400,  500,  and  600  psia. 
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FIGURE  23.  Variation  of  the  measured  transition  Reynolds  number  with 
the  unit  Reynolds  number. 


FIGURE  26.  Test  set-up  for  the  wake  turbulence  experiment. 
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THEORY  (REF.  2) 


FIGURE  30.  Minimum  stagnation  temperature  required  for  possible 
FSL  flows  if  and  \  are  first  specified. 


APPENDIX  A 


THE  DIVIDING- STREAMLINE  TEMPERATURE  IN  SHEAR  LAYERS 
ACCORDING  TO  THE  CROCCO  RELATION 

The  objective  of  this  Appendix  is  to  derive  an  expression  for  the  static 
temperature  TDSL  on  the  dividing  streamline ,  which  accounts  for  Mach  number 
and  speed  ratio  variation  from  low  to  high  levels,  as  well  as  arbitrary  total 
temperature  differences  between  the  two  streams.  Excepting  differences  in 
molecular  composition,  this  would  then  provide  the  information  needed  to  en¬ 
compass  all  conceivable  physical  ranges  of  free  shear  layers. 

Our  task  is,  actually,  to  simplify  and  exploit  the  temperature-velocity 
relation  ("Crocco  relation")  derived  by  Shapiro  (Ref.  24),  Korst  and  Chow 
(Ref.  25)  and  others: 
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(A.  1) 


The  objective  here  is  to  evaluate  this  formula  at  the  dividing  streamline 
(or  "DSL"): 
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U1 


A  U- 


(A.  2) 


and  express  the  result  in  terms  of  M^  X  and  T02/Tq1.  Note  that 

f(Mt)  -  Mi2 

To  do  this,  note  that 
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These  equations  are  still  unsatisfactory  because  they  do  not  contain  the  ratio 
t02/T0i  explicitly;  further  manipulation  produces  the  desired  result: 


..  .  V-l  „  2.  /1+l02.  (Y  "  X)  M1 

T,  2  (1  2  M1  )  (  T01}  2(1+X)  2 


_DSL  .  1 

ri 


(A.  9) 


In  equation  (A. 9)  the  parameters  X  and  Tq  /Tq1  are  independent,  and  any 
set  of  them  should  define  a  unique  T^j/Tj  value.  Physically  speaking,  how¬ 
ever  t  there  is  an  exception  made  when  the  second  term  on  the  r.h. s.  of  (A. 9) 
exceeds  the  first  which,  it  will  be  noted,  can  happen  especially  if 
T02/toi  <  <  1*  In  this  case  TDgL  <  0,  a  clearly  impossible  (but  algebraical] 

permitted)  result.  There  are,  therefore,  minimum  allowed  values  of  Tr„/T, 


permitted)  result.  There  are,  tnererore,  minimum  anowea  vawes  *02  01 

which,  combined  with  any  pair  of  prescribed  ML, A  values,  will  give  T  >  0 
according  to  (A. 9),  and  these  are  plotted  in  Figure  30.  An  additional,  but 
similar  restriction  on  Tqo/Iqi  arlses  when  ^  re<lulre  that  in  the  yPothe" 
sized  ''slow”  stream  (subscript  "2")  the  temperature  T-  >  0.  This  criterion 
can  be  derived  from  Eq.  (A. 7)  and  (A. 8)  and  is  plotted  on  the  same  figure. 

The  criterion  required  to  keep  Tg  >  0  is,  as  Figure  30  shows,  more  restrictive 
than  the  requiring  TDSL  >  °»  tha^  is,  the  latter  criterion  can  be  ignored  so 

long  as  the  former  is  accounted  for.  Thus,  according  to  Figure  30,  if  we 
hypothesize  a  FSL  with  M,  -  3  and  A-  0.1,  say,  then  we  can  prescribe  to  this 
flow  any  T02/Tq1  larger  than  0.42  or  so.  Values  of  TQ2/T0l  less  than  that 

will  create  a  physically  unattainable  flow.  Thus,  Equation  (A. 9)  subject  to 
the  restrictions  of  Figure  30  is  the  necessary  tool  for  computing  the  kine¬ 
matic  viscosity  entering  the  transition  formula  in  the  text. 

We  will  forego  demonstrating  the  lateral  variation  of  temperature  across 
the  FSL  according  to  (A. 1),  since  such  plots  are  relatively  easy  to  do.  It 
is,  however,  important  to  keep  in  mind  that  there  are  values  of  Tq2' Tqj  I°r 


which  T  does  not  vary  monotonically  from  T^  to  T2< 
that 


This  can  be  seen  by  noting 
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and  plotting  (T-T^/Tj  va. 


T^where.  as  usual: 
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The  result  will  show  that,  especially  for  large  M.  and  X  and  for  <  1» 

T  arrives  at  a  —»<■■»  between  the  two  streams.  The  physical  reason,  of  course, 
is  that  as  the  fast  stream  slows  down  to  match  the  speed  of  the  slow  one,  it 
experiences  an  isentroplc  temperature  rise.  In  such  cases  the  temperature  pro¬ 
file  resembles,  as  it  should,  a  cold-wall  supersonic  boundary  layer. 


This  parenthetic  remark  was  made  to  prepare  the  user  of  Eq.  (A. 9)  for  oc¬ 
casional  "odd"  behavior  of  TDg, ,  such  as  an  actual  decrease  for  certain  con¬ 
ditions  illustrated  on  Figure  31  on  which  Eq.  (A. 9)  is  plotted.  Such  behavior, 
setting  in  for  T02/T0-  <  1,  will  have  consequences  in  the  transition  predictions 
found  in  that  ten. 


